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Abstract 

Travel  time  perturbations  of  adiabatic  normal  modes  due  to  an  internal  tide  and 
internal  mode  field  in  the  Barents  Sea  are  examined.  A  formalism  for  the  travel  time 
perturbation  due  to  a  change  in  sound  speed  is  presented.  Internal  tide  and  internal 
wave  amplitude  spectra  are  calculated  from  Brancker  temperature  loggers  which  were 
deployed  on  moorings  in  the  Barents  Sea  during  the  August  1992  Barents  Sea  Polar 
Front  Experiment.  In  particular,  the  first  three  internal  wave  mode  amplitudes  are 
estimated  from  the  four  Brancker  temperature  loggers  on  the  southwest  mooring  of 
the  array.  Modal  perturbations  in  acoustic  pulse  travel  time  and  the  travel  time 
covariance  are  calculated  and  compared  for  consistency  to  a  simple  ray  model.  These 
perturbations  are  small  for  the  modal  arrivals  that  the  vertical  acoustic  array  which 
was  deployed  is  expected  to  resolve.  The  third  internal  wave  mode  has  the  largest 
impact  on  the  acoustic  arrivals,  per  unit  amplitude,  but  the  first  internal  wave  mode 
dominates  the  scattering  due  to  having  a  much  larger  amplitude  overall. 
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Chapter  1 


Introduction 


Tomography  comes  from  the  Greek  word  tomos  meaning  slice  or  section.  Ocean 
acoustic  tomography  is  the  science  of  taking  acoustic  slices  or  sections  of  the  ocean 
sound  speed,  temperature,  and  current  fields.  In  this  thesis  we  will  deal  with  the 
application  of  tomographic  techniques  to  the  study  of  internal  waves  and  internal 
tides. 

1.1      Motivation 

The  ocean  covers  over  two-thirds  of  the  earth's  surface.  With  its  tremendous  heat 
capacity  and  enormous  currents,  the  ocean  exerts  a  significant  influence  on  the  world's 
climate  and  meteorology.  This  influence  is  not  well  understood,  in  part  due  to  the 
difficulty  in  observing  the  ocean  volume,  and  also  in  part  due  to  the  ocean's  meso- 
scale  and  inter-annual  variability.  Remote  sensing  using  electromagnetic  energy  in 
the  forms  of  light  and  radar,  which  has  allowed  man  to  peer  billions  of  light-years 
into  outer  space,  encounters  in  the  ocean  an  opaque  medium.  The  average  depth  of 
the  ocean  is  4000  m,  but  blue-green  light  and  radar  can  only  penetrate  to  depths 
of  100  m  and  the  surface,  respectively.  What  knowledge  we  have  about  the  ocean 
depths  was  obtained  mainly  through  in-situ  measurements  conducted  by  a  myriad 
of  research  vessels  performing  tens  of  thousands  of  casts  of  Nansen  bottles,  Niskin 
bottles,  and  CTD's. 
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Though  bottle  and  CTD  casts  are  important  and  still  continue  today,  they  are 
time  and  labor  intensive.  Not  only  are  they  costly,  but  the  slow  sampling  rate  makes 
them  susceptible  to  aliasing  of  many  oceanic  processes.  As  new  methods  have  been 
sought  out,  the  use  of  acoustics  to  probe  the  depths  of  the  ocean  has  become  an 
attractive,  alternate  technology  in  oceanic  research.  Acoustics  has  allowed  man  to 
determine  accurately  the  topography  of  the  ocean  floor,  to  track  marine  life,  to  profile 
the  ocean  currents,  and  to  perform  many  other  useful  tasks.  One  of  the  more  recent 
research  efforts  in  acoustics  has  been  that  of  ocean  acoustic  tomography.  The  idea 
behind  tomography  is  to  take  many  slices  or  sections  of  an  object  at  different  angles  to 
produce  a  three  dimensional  picture  of  the  subject.  The  theory  was  mathematically 
proven  by  Radon  in  1917  [cited  in  Munk  and  Wunsch,  1979].  Tomography  has  been 
widely  used  in  medicine  (CAT  scan  stands  for  Computer  Assisted  Tomography  scan), 
and  in  seismology. 

The  basic  principle  behind  ocean  acoustic  tomography  is  that  by  measuring  the 
long-term,  temporal  variation  in  acoustic  pulse  travel  times  between  a  source  and 
a  receiver,  due  to  the  sound  speed  field  variation  from  effects  such  as  fronts  and 
eddies,  one  can  image  these  features  acoustically,  and  thus  study  them.  This  temporal 
variation  had  been  noted  in  several  early  experiments  [see  Hamilton,  1977;  Steinberg 
and  Birdsall,  1966  (both  cited  in  Munk  and  Wunsch,  1979)],  but  was  often  considered 
noise.  Worcester  [1977]  conducted  one  of  the  first  experiments  in  ocean  acoustic 
tomography,  when  he  measured  reciprocal  acoustic  travel  times  between  two  ships  25 
km  apart,  in  order  to  investigate  the  feasibility  of  using  the  travel  time  differences 
to  monitor  mid-ocean  currents.  Munk  and  Wunsch  [1979]  proposed  expanding  the 
range  between  source  and  receiver  to  1000  km  to  monitor  mesoscale  ocean  features. 
Since  then,  acoustic  tomography  has  been  applied  to  both  small  scale  and  large  scale 
oceanographic  eddies,  as  well  as  to  shallow  water  and  deep  water  environments. 
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1.2  Thesis  Objectives 

In  this  work,  we  will  first  calculate  internal  tide  and  internal  wave  spectra  and  am- 
plitudes from  temperature  measurements  taken  at  moorings  at  several  depths  in  the 
Barents  Sea  during  the  August  1992  Barents  Sea  Polar  Front  (BSPF)  Experiment. 
Using  this  oceanographic  data,  the  perturbation  in  acoustic  travel  times  caused  by 
the  internal  tides  and  internal  waves  is  predicted.  This  will  be  eventually  compared 
to  actual  perturbations  in  travel  time  recorded  during  the  experiment.  These  pertur- 
bations are  expected  to  be  significant  enough  to  be  measurable. 

1.3  Thesis  Content 

The  scope  of  this  thesis  is  as  follows:  Chapter  two  first  presents  the  theory  behind 
acoustic  tomography,  internal  tides  and  internal  waves.  Next,  using  adiabatic  normal 
mode  theory,  a  formulation  for  the  travel  time  perturbation  of  an  acoustic  signal  due 
to  internal  waves  is  presented.  Chapter  three  describes  the  region  of  the  experiment, 
the  experiment  itself,  and  the  data  collected.  Chapter  four  details  the  oceanographic 
data  analysis  conducted  and  reviews  the  assumptions  which  have  been  made.  Chap- 
ter five  examines  travel  time  perturbations  due  to  internal  tides  and  internal  waves. 
Comparisons  between  a  simple  ray  model  and  the  modal  formulation  derived  in  Chap- 
ter two  are  made.  Finally,  Chapter  six  presents  conclusions  and  gives  suggestions  for 
future  work. 
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Chapter  2 


Theory 


2.1      Acoustic  Tomography 

Acoustic  tomography  can  be  broken  into  a  forward  and  an  inverse  problem.  In  the 
forward  problem,  the  travel  time  of  a  pulse  along  a  path  is  computed  from  the  sound 
speed  and  current  fields  of  the  ocean.  In  the  inverse  problem,  given  a  set  of  travel 
times,  the  sound  speed  and  current  structure  of  the  ocean  is  to  be  estimated.  In  this 
section,  we  will  look  at  the  ray  theory  approach  to  the  forward  problem. 

The  formulation  of  the  forward  problem  given  here  follows  that  of  Howe  et  al. 
[1987].  Using  geometric  optics,  the  travel  time  T,  along  a  ray  path  T,  in  an  ocean 
having  a  sound  speed  field  C(x,t)  and  a  current  field  u(x,tf)  is 

T  =  f  - (2  1) 

'      Jr,  C(x,0  +  u(x,<)-t'  V  '  ; 

where  s  is  the  arc  length,  and  r  is  a  unit  vector  tangent  to  the  ray.  The  sound  speed 
field  C(x,<)  can  be  decomposed  as 

C(x,t)  =  C0(x)  +  £C(x,t),  (2.2) 

where  Co(x)  is  a  "known"  background  sound  speed  field  and  6C(x,  t)  is  a  sound  speed 
perturbation  such  that  8C(x,t)  <C  Co(x).  From  these  equations  it  can  be  seen  that 
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a  ray  traveling  with  a  current  has  a  smaller  travel  time  than  a  ray  traveling  against 
the  current. 

The  perturbation  travel  times  for  a  reciprocal  transmission  between  transceivers 
can  be  found  by  substituting  (2.2)  into  (2.1)  and  linearizing.  This  gives 

ST?  =  n-T«  =  - ^  ^^c^^^  (2-3) 


and 

STt    =  Tx    -  To,  =  -  J^ ^- da,  (2.4) 

where  the  +  superscript  refers  to  a  ray  propagating  in  the  -f  x  direction  and  the 

—  superscript  is  a  ray  propagating  in  the  — x  direction.    By  taking  the  sum  s,,  or 

the  difference  d{  of  (2.3)  and  (2.4),  either  u(x,£)  •  r,  or  6C(x,t)  may  be  eliminated, 

respectively  giving 

26C(x,t) 


s,  =  ST?  +  ST-  =  -fr^ds,  (2.5) 


and 


* 


-^-^=-£w*-  (26) 


Howe  et  ah  [1987]  showed  that  internal  waves  and  large  scale  velocity  shears  do 
not  significantly  perturb  the  reference  ray  path  so  that  for  a  two-way  transmission  the 
ray  paths  are  reciprocal  to  first  order  and  for  a  one-way  transmission  the  ray  paths  for 
the  reference  ray  and  the  perturbed  ray  are  identical.  In  the  BSPF  experiment,  only 
the  one-way  transmission  case  is  examined.  By  noting  that  6C(x,t)  >>  u(x,£)  •  r, 
(2.3)  can  be  reduced  to 

ST?  =  Tt-T0>  =  -Jr^ds.  (2.7) 


15 


2.2      Internal  Motions  in  the  Ocean 

Internal  motions  in  the  ocean  of  interest  in  this  thesis  are  internal  tides  and  internal 
waves.  This  classification  is  a  matter  of  the  frequency  of  interest.  Internal  tides  are 
forced  at  tidal  frequencies  u>tide,  and  are  either  propagating  or  evanescent  depending 
on  whether  Utidt  is  greater  than  or  less  than  the  inertial  frequency  /,  respectively.  The 
propagating  internal  wave  band  is  between  the  inertial  frequency  /  and  the  highest 
buoyancy  frequency  Nmax.  These  internal  motions  are  clearly  seen  in  Fig.  2-1,  which 
is  an  example  of  a  temperature  record  that  has  been  has  been  Fourier  transformed 
into  its  frequency  components  (using  methods  to  be  discussed  later  in  Chapter  4). 
Clearly  shown  in  the  figure  are  the  semidiurnal  internal  tide  along  with  its  harmonics, 
and  the  internal  wave  spectrum  at  higher  frequencies.  Nmax  at  this  depth  is  215  cycles 
per  day. 

2.2.1      Internal  Tides 

Internal  tides  are  generated  by  the  surface  or  barotropic  tides  advecting  water  over 
the  ocean  bottom.  Any  change  in  the  bottom  topography  causes  the  isopycnals  to  be 
displaced.  This  disturbance  excites  the  internal  modes  in  the  water.  The  continental 
shelf  break  is  a  principal  generation  site  of  internal  tides,  whose  generation  mechanics 
have  been  well  studied  [e.g.  Baines,  1974;  Rattray  et  al.  1969].  Internal  tides  can 
also  be  generated  at  other  topographic  changes  such  as  guyots,  troughs,  sills,  and  sea 
mounts  [e.g.  Maxworthy,  1979]. 

Wunsch  [1975]  reviewed  several  internal  tide  records  and  noted  they  often  showed 
a  quasi-semidiurnal  or  quasi-diurnal  periodicity,  but  that  due  to  the  record  length, 
background  noise,  and  frequency  shifting  of  the  tidal  peak,  it  was  often  hard  to  resolve 
the  tidal  constituents.  He  also  noted  that  the  tidal  signal  was  intermittent.  Baines 
[1986]  suggested  that  the  intermittent  nature  of  the  internal  tide  might  be  due  to 
the  changes  in  stratification  between  the  generation  site  and  the  observer,  or  due  to 
frontal  changes  or  upwelling.  The  generation  of  internal  tides  is  a  very  non-linear 
effect.    Many  records  show  periodicities  of  6  hours  (twice  semidiurnal)  or  8  hours 
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Figure  2-1:  Power  Spectrum  Density  of  Brancker  temperature  logger  #3663  located 
on  NE  array  at  a  depth  of  32.1m.  Evident  are  the  M2  tide,  its  higher  harmonics,  and 
the  internal  wave  band.  Nmax  at  this  depth  is  215  cpd.  The  estimate  has  4  degrees 
of  freedom. 

(diurnal  plus  semidiurnal),  or  even  3  hours  [Baines,  1986].  The  3  hour  (8  cpd),  6 
hour  (4  cpd),  and  12  hour  (2  cpd)  periods  are  evident  in  Fig.  (2-1).  Some  of  these 
higher  harmonics  of  M2  could  be  due  to  computing  the  spectrum  of  temperature, 
rather  than  density,  however,  as  will  be  discussed. 

In  his  review,  Wunsch  [1975]  also  noted  that  most  of  the  internal  tidal  energy  is 
contained  in  the  first  few  lowest  modes  with  mode  one  dominant.  The  direction  of 
propagation  is  generally  shoreward,  but  this  need  not  be  the  case.  For  instance,  Levine 
and  Richman  [1989]  performed  a  study  which  showed  that  the  internal  tide  off  the 
coast  of  California  propagated  along  the  shelf.  They  suggested  that  the  internal  tides 
were  either  being  refracted  from  fronts  along  the  coast  or  from  mesoscale  structures, 
or  that  the  resultant  direction  is  due  to  many  waves  being  radiated  from  several  shelf 
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breaks  oriented  in  different  directions.  Another  study  by  Holloway  [1984]  found  that 
the  first  mode  dominates  and  propagates  onshore  at  an  angle  of  about  30  °  from 
normal  to  the  bathymetry. 

2.2.2      Internal  Waves 

Internal  waves  are  found  throughout  the  ocean.  They  exist  in  the  ocean  interior 
due  to  disturbances  between  fluids  of  differing  densities.  Although  numerous  internal 
wave  generation  mechanisms  have  been  detailed,  only  the  more  important  mechanisms 
will  be  discussed  here.  There  exist  several  thorough  reviews  which  the  reader  may 
wish  to  consult  [Thorpe,  1975;  Muller  and  Olbers,  1975;  LeBlond  and  Mysak,  1978; 
Olbers,  1983].  The  internal  wave  generation  mechanisms  involve  interactions  with  the 
bottom,  the  surface,  or  the  interior. 

Bottom  topography  is  an  important  mechanism  in  the  formation  of  internal  waves, 
particularly  in  coastal  regions.  This  mechanism  has  been  extensively  studied  in  the 
region  of  the  continental  shelf  break  [Sawyer,  1983;  Curtin  and  Mooers,  1975],  over 
ridges  in  bays  [Halpern,  1971;  Chereskin,  1983]  and  over  sills  in  straits  [Gargett,  1976; 
Farmer  and  Smith,  1980].  In  a  horizontally  stratified  ocean,  the  flow  of  the  baroclinic 
tide  over  a  sharp  boundary  (such  as  a  continental  shelf  or  a  ridge)  causes  a  shear 
instability  and  creates  a  lee  wave.  Depending  upon  the  specific  conditions,  the  lee 
wave  will  either  propagate  shoreward  as  a  lee  wave  train  or  evolve  into  an  undular 
bore  structure.  These  internal  wave  packets  or  solitons  have  generally  been  observed 
when  the  tide  turns. 

Surface  interactions  can  also  generate  internal  waves.  Thorpe  [1975]  reviewed 
several  surface  interaction  mechanisms,  which  we  present  here  in  order  of  assumed 
importance.  First,  the  wind  stress  can  produce  a  traveling  stress  field  in  the  mixed 
layer,  which  results  in  a  divergent  current  field  in  the  mixed  layer.  The  vertical 
Ekman  velocities  created  by  the  wind  stress  can  excite  internal  waves.  Second,  a 
traveling  pressure  field,  created  by  atmospheric  pressure,  can  resonantly  couple  with 
the  internal  waves  to  form  waves  that  grow  linearly.  Third,  a  traveling  buoyancy  flux, 
produced  by  solar  heating  and  precipitation,  induces  a  vertical  velocity  at  the  lower 
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boundary  of  the  surface  layer.  The  induced  velocities  excite  internal  waves.  Fourth, 
a  spectrum  of  surface  gravity  waves  can  be  produced  by  the  wind  stress  acting  over 
the  ocean  surface.  Two  of  these  higher  frequency  surface  gravity  waves,  si  and  s2 
can  beat  against  each  other  to  produce  a  lower  frequency  internal  wave  with  wave 
number  ktu, 

kiw  —  ksi  —  kS2,  (2-8) 

and  frequency  Uiw 

Uiw  =  w,i  -  ua2.  (2.9) 

The  interaction  of  the  surface  waves  generates  high  frequency  internal  waves  in  the 
upper  ocean,  but  this  mechanism  is  insignificant  for  the  deep  ocean.  Under  conditions 
of  a  rough  sea  and  a  strong,  shallow  thermocline,  however,  it  may  be  comparable  to 
the  wind  stress  mechanism  (the  first  mechanism  mentioned)  [Olbers,  1983].  Finally, 
instability  between  the  surface  boundary  layer  and  the  deep  ocean,  such  as  turbulence 
in  the  mixed  layer,  may  also  generate  internal  waves. 

As  regards  the  role  of  the  interior  of  the  ocean  in  the  generation  of  internal  waves, 
very  little  is  known.  LeBlond  and  Mysak  [1978]  list  as  possible  generation  mechanisms 
processes  such  as  penetrative  convection,  geostrophic  adjustment,  and  internal  tidal 
modulation  of  the  thickness  of  the  thermocline. 

Garrett  and  Munk  [1972]  looked  at  several  internal  wave  observations  and  created 
from  them  a  model  of  the  internal  wave  field  frequency- wavenumber  spectrum.  The 
Garrett- Munk  model  assumes  horizontal  isotropy  and  vertical  symmetry.  It  has  the 
form 

S(kh,u>)  oc  ^-^-p+1(^2  -  J2)-1'2  (2.10) 

where  kh  is  the  horizontal  wavenumber,  /  is  the  local  inertial  frequency,  and  fi  is  the 
cut-off  wavenumber  which  satisfies 

jiocw'-1^2-/9)1'8.  (2.11) 

Substituting  in  (2.11)  into  (2.10),  assuming  a,'2  >  /2,  and  then  integrating  (2.10)  over 
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kh,  gives  for  the  frequency  spectra 

S(u)  =  u~p.  (2.12) 

The  indices  p  and  r  are  determined  from  observed  spectra,  which  suggested  values 
of  p  between  5/3  and  2,  and  r  between  2/3  and  3/2.  Garrett  and  Munk  proposed  a 
value  of  2  for  both  r  and  p.  The  Garrett-Munk  model  has  been  modified  and  refined 
over  the  years  as  new  experimental  information  has  become  available. 

While  the  spectrum  for  the  deep  ocean  is  well  modeled  by  Garrett-Munk,  that 
for  shallow  waters  is  not.  An  internal  wave  experiment,  conducted  as  part  of  the 
GARP  Atlantic  Tropical  experiment  or  GATE  in  1974,  gave  a  spectrum  which  devi- 
ated strongly  from  the  Garrett-Munk  model  [Kase  and  Siedler,  1980].  The  frequency 
spectrum  showed  three  dominant  energy  bands  centered  about  the  near-inertial  fre- 
quency, the  M2  tidal  frequency,  and  3  cycles  per  hour  (cph).  The  high  frequency 
wave  field  was  intermittent  and  anisotropic.  The  non-stationarity  of  the  spectra  is  a 
result  of  the  variability  of  the  atmospheric  generation  mechanisms  [Kase  and  Clarke, 
1978].  Kase  and  Clarke  also  argue  that  the  the  higher  frequency  waves  are  trapped 
near  the  surface  where  a  maximum  A^  overlies  a  thick  weakly  stratified  lower  layer. 
These  waves  form  a  low  order  standing  mode.  Baines  [1986]  also  noted  a  study  in 
which  for  depths  less  800m,  25-45%  of  the  energy  was  in  mode  1  which  propagated 
shoreward. 

Dynamical  Equations  for  Internal  Waves 

Using  the  basic  equations  of  momentum  and  continuity,  an  exact  modal  solution  for 
internal  waves  can  be  derived.  The  equation  of  momentum  for  an  isentropic  fluid  in 
the  absence  of  viscous  effects  is  given  by  [e.g.  Gill,  1982] 

Du 

—  +  2nxu  =  -/T1Vp-g.  (2.13) 
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and  the  continuity  equation  for  an  incompressible  fluid  is  given  by 

V  -u  =  du/dx  +  dv/dy  +  dw/dz  =  0.  (2.14) 

Following  the  derivation  of  Apel  [1987],  a  small  perturbation  in  p  and  p  is  applied  to 
(2.13).  The  perturbation  pressure  and  density  are  then  given  respectively  by 

p=-gpz  +  p',  (2.15) 

and  by 

p  =  Po  +  p'.  (2.16) 

The  momentum  equation  (2.13)  is  linearized  by  neglecting  the  products  of  perturba- 
tion quantities,  and  becomes 

du       „  1  dp' 

m-fv—7o^'  (2'17) 

!-/«  =  -M^  (2.18) 

dt  p0  dy 

dw  1  dp'       p' 

at         p0  dz      pQ 

where  /  is  the  inertial  or  "local  Coriolis"  frequency.    We  now  assume  a  separable 
solution  for  the  vertical  velocity  w(x)  of  the  form 

io(x,  t)  =  w(z)  exp[i{kx  +  ly  -  tot)},  (2.20) 

where  w(z)  is  the  vertical  eigenfunction.     The  vertical  velocity  u>(x)  satisfies  the 
boundary  conditions  of  no  normal  flow  at  the  bottom 

w(x)  =  0    at  z  =  h  (bottom),  (2.21) 
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and  that  a  particle  on  the  free  surface  remains  on  it 


w(x)  =  drj/dt    at  z  =  tj  (surface), 


(2.22) 


which  can  be  applied  at  z  =  0  since  77  is  small.  This  approximation  introduces  little 
error,  as  the  vertical  disturbance  at  the  air-water  boundary  is  much  smaller  than  the 
maximum  amplitude  by  an  order  of  1000.    By  manipulating  equations  (2.14-2.19), 

one  can  obtain 

'N2(z)-u2' 


d2w 


dz< 


+  « 


u;2-/2 


w  =  0, 


(2.23) 


where  kk  is  the  horizontal  wavenumber  that  satisfies 


k2  =  k2  +  /2, 


(2.24) 


and  where  N  is  the  buoyancy  or  Brunt- Va.isa.la  frequency  which  satisfies 


N2  =  - 


9  dp0 
Po  dz 


(2.25) 


For  u>  <  f  or  u>  >  N,  the  vertical  wave  number  is  imaginary  and  the  solutions  to 
(2.20)  are  evanescent. 

By  using  the  relationship 


w  = 


dV 


where  C,  is  the  particle  displacement,  and  substituting  in  (2.23)  gives 


d2(       2 

dz~2  +  kh 


N2{z)-u< 


u2-P 


C  =  o, 


with  the  following  boundary  conditions: 


(2.26) 


(2.27) 


C,{z)  =  0    at  z  =  0  (surface); 


(2.28) 


and 


£(z)  =  0    at  z  =  h  (bottom). 
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(2.29) 


This  equation  can  be  formulated  as  an  eigenvalue  equation  and  is  easily  solved  using 
numerical  methods  such  as  finite  element  methods.  The  exact  modal  solution  is 

oo 

(  =  Yi  MM  exp  [i  (k,x  +  l,y  -  ut)} ,  (2.30) 

where  /  =  1,2,...  is  the  mode  number,  Ai  is  the  amplitude,  and  tpi(z)  is  the  vertical 
internal  wave  mode  eigenfunction. 

2.3     Adiabatic  Acoustic  Normal  Modes 

We  now  return  to  acoustics.  In  a  range  independent  environment,  the  acoustic  pres- 
sure p(z,x)  for  a  point  source  at  r  =  0,  z  =  Z\  is  described  by  the  Helmholtz  equation 

V2p  +  K\z)V  =  —6{z  -  z1)S{r).  (2.31) 

r 

Brekhovskikh  and  Lysanov  [1991]  give  as  the  solution 

P(r,z)  =  TTlJ^trnizWmiz^H^iKrnr).  (2.32) 

m 

The  vertical  mode  function  <pm  satisfies 


dz 


a  + 


K\z)  -  A'£   4>n  =  0,  (2.33) 


where  K{z)  is  the  total  wavenumber  and  A'^  is  the  mih  modal  eigenvalue.  Equa- 
tion (2.32)  is  a  good  approximation  for  many  parts  of  the  ocean,  but  as  it  is  range 
independent,  it  is  inadequate  where  the  depth  of  the  ocean  changes  significantly, 
where  the  sound  waves  cross  ocean  fronts,  or  where  there  is  a  significant  change  in 
the  sound  speed  profile  along  the  propagation  path.  In  these  cases  the  Helmholtz 
equation,  equation  (2.31),  must  be  modified  to  take  into  account  the  variation  of 
ocean  parameters  with  range  and  becomes 

V2p  +  K2{z,r)P  =  Q.  (2.34) 
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Shang  [1989]  gives  as  the  solution  to  (2.34) 

p(r,z)  =  f^)^),  (2.35) 

n=l  Vr 

where  Fn(r)  is  the  radial  component  and  the  vertical  mode  function  <f>n(z,r)  satisfies 

dz  \po{z)  dz  J       \    p0(z)  p0{z)  J 

The  mode  functions  (f>n  are  mutually  orthogonal  and  satisfy 


fh   l 

/      —(f>m(f>n  =  £>mn-  (2.37) 

JO     Pn 


Substituting  (2.35)  into  (2.34),  multiplying  by  <f>m,  and  then  integrating  over  z  from 
0  to  h  results  in  the  following  coupled  equation 

^f  +  (*iW  +  ^-)Fm  =  -E  [Amn(r)Fn  +  2Bmn(r))  ^,  (2.38) 

where  Amn  and  Bmn  are  coupling  coefficients  defined  as 

Amn  =  Jp^<f>m(r,z)^(z,r)dz,  (2.39) 

Bmn  =  j  PolK(r,z)^-(z,r)dz.  (2.40) 

To  obtain  the  "adiabatic  approximation",  we  neglect  the  l/4r2  term  in  (2.38)  since 
we  are  interested  in  solutions  with  Kmr  ^>  1,  and  we  assume  that  the  mode  coupling 
effects  are  small.  Physically,  this  means  that  the  ocean  properties  (depth,  sound 
speed)  "vary  slowly"  as  a  function  of  horizontal  range.  The  phrase  "vary  slowly" 
means  that  the  fractional  change  of  the  ocean  properties  over  a  distance  of  less  than 
a  wavelength  is  small  [Pierce,  1965].  This  type  of  medium  has  been  termed  "almost 
stratified".  The  adiabatic  approximation  solution  to  (2.34)  is 

p(r,z)  =  v^exp(?7r/4)  V^%J£lM^exp  (i  [*  Km(r')dA  .  (2.41) 

m         y/Km(r)r  V  Jo  J 
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In  the  Barents  Sea,  the  background  sound  speed  profile  and  the  depth  vary  rea- 
sonably slowly  along  the  sound  propagation  path  on  either  side  of  the  front.  Taking 
advantage  of  this  we  can  crudely  treat  the  path  on  either  side  of  the  front  as  a  "nearly 
range  independent"  segment,  and  use  the  adiabatic  approximation.  We  expect  strong 
coupling  of  the  modes,  though,  in  going  through  the  frontal  region. 

2.4      Sound  speed  perturbations   due  to  Internal 
Waves 

Equation  (2.3)  shows  that  the  two  mechanisms  affecting  the  perturbation  of  the  acous- 
tic signal  are  the  horizontal  current  u(x,  t)  •  r  and  the  sound  speed  perturbation 
6C(x,t)  which  is  related  to  the  particle  displacement  £  using  the  following  relation- 
ship 

SC  =  ^(,  (2.42) 

oz 

where  C'p  =  dCv/dz  is  the  potential  sound  speed  gradient  and  (  is  the  solution  to 
equation  (2.27).  Assuming  that  both  mechanisms  are  due  to  internal  waves,  Essen  et 
al.  [1983]  showed  that  SC  has  a  much  greater  affect  on  the  acoustic  signal  than  does 
u(x,<)  •  r. 

Following  a  recent  derivation  by  Jin  [1992],  for  a  band  of  frequency  u;,  equa- 
tion (2.30)  can  be  generalized  to 

£  =  Xw      /     ai{u,0)i)l{u,z)exv[i{Hu)x  +  ll{u)y -ut)](L}dO,  (2.43) 

where  <n(u,6)  is  the  internal  wave  amplitude  spectrum  for  mode  /,  xp[(u,z)  is  the 
internal  wave  mode  function,  0  is  the  azimuthal  angle,  and  ki(u)  and  h(u))  are  the  in- 
ternal wave  eigenvalues  that  satisfy  equation  (2.24).  The  quantities  iI)i(uj,z),  a/(u>,#), 
ki(to),  and  //(u>)  are  assumed  to  be  range  independent. 

We  next  follow  Rajan  et  al.  [1987]  and  perturb  the  background  acoustic  wavenum- 
ber  A'°(r,  z)  by  perturbing  the  sound  speed  profile.  A  new  wavenumber  A'(r,  z)  can 
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then  be  written  as 


K{r,z)  =  u>/[C0(r,z)  +  6C{r,z)]. 


(2.44) 


Using  standard  perturbation  theory,  a  perturbation  eigenvalue  6Km(r)  for  (2.41)  is 
found  to  be 


6K^  =  ^77Wr-J      ^)(^z)\SK(r,z)2\<f>^(r,z)dz.  (2.45) 


2K£'{r) 


The  middle  factor  in  the  integral  in  (2.45)  can  be  rewritten  as 


6K{r,z)2    =     K2{r,z)-K(0)2{r,z) 


ur 


C$(r,z) 


-2 


rW3)  - 


S     -Km{r,z){26C/Co) 


(2.46) 


Substituting  (2.46)  into  (2.45)  results  in 


6Km(r)  =  -±-  /V'M  U&(,,x)f  (*<0>(r,*))'  ^^<i2.  (2.47) 

/<V(r)    °  i    v  '     C0(r,2) 


Using  the  WKB  approximation,  the  adiabatic  modal  phase  is 


(p{r)  =  /     Km(r)dr, 
Jo 


(2.48) 


and  the  background  modal  phase  is 


<f°(r)  =  [R  K°m(r)dr  =  K°m 
Jo 


(2.49) 


The  modal  phase  perturbation  is  thus 


fR 
<Vm  =  /     6Km(r)di 
Jo 


(2.50) 
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If  the  received  signal  is  processed  by  a  narrow  band  filter  with  center  frequency 
u>o,  then  the  group  travel  time  of  each  modal  pulse  is  given  by  the  stationary  time  as 


'm    — 


d[ip, 


du> 


(2.51) 


so  that  the  group  travel  time  perturbation  is 


6tm  = 


d[8tp 


m]u>o 


<JLo 


(2.52) 


Inserting  (2.47)  into  (2.50),  and  (2.50)  into  (2.52)  and  using  the  results  of  equa- 
tion (18)  of  Lynch  et  al.  [1991]  results  in 


i-R     roo 

Stm=         /     Gm{r,z)6C(r,z)drdz, 
Jo    Jo 


(2.53) 


where  Gm(r,  2),  the  group  velocity  dispersion,  is 


Gm{r,z)  =  p0 


-1 


1^(0)^3 
Jim    ^0    L 


Ct0)M 


r  ~  w$* 


)|A°>(r,z)f  +  u,*A|40)(r)Z) 


,    (2.54) 


and  Cj°\r)  and  Uffl(r)  are  the  modal  phase  and  group  velocities,  respectively,  for 
the  unperturbed  case.  By  assuming  range  independence,  and  also  that 


uL°\r)  ~    ' 


(2.55) 


(2.54)  becomes 


Gm{z)  -  -p^ 


1 
ni 


4&(z) 


+  Mo- 
dus 


e\z) 


(2.56) 


where  p(z)  =  1  and  C$  «  Co-  The  bar  denotes  the  average  sound  speed. 
The  covariance  of  Stm  is  formed  by 


covm{R,  t)  =  (6tm{R,t)8t*m{R,  t  +  r)), 


(2.57) 
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which  using  (2.53)  becomes 


-R    rR 
cov 


Jrri     rtt     roo     roo 
'      /      /      /      Gm(r,z)Gm(r',z')x 
0    Jo    Jo     Jo 

(6C{r,  z,  t)6C{r',  z',  t  +  r))drdr'dzdz\  (2.58) 


where  from  (2.42)  and  (2.43) 


00      rN     r2n 


6C(x,y,z,t)    =    C'p(z)Y,         /      a,(u>,0)1,l(u>,z)x 

jr{  J i  Jo 

exp[i(khi{u)x  cos(O)  +  k  hi  {u)y  sin(0)  —ujt]dujd9.        (2.59) 

Assuming  that  the  amplitudes  a/(u>,0)  are  independent,  random  variables  with  zero 
means  implies  that 

(ai(u,0))  =  O, 

and 

(ai(u,0)ai{w,O))  =  \(h{lo,6)\2  6ww,6ee>6u>, 

where  6^1,  8gg>,  and  <$///  are  three  delta  functions.  Additionally,  if  it  is  assumed  that 
the  x  axis  is  in  the  direction  from  the  source  to  receiver,  we  obtain 

{8C(r,z,t)6C(r\z',t  +  T))  = 

°°       rN     r2ir 

jr[Js  Jo 

C'p(z)C'p{z')  exp  [i{kM(u>)  cos  0{r-  r')  -  wr)]  diodd.  (2.60) 

Substituting  (2.60)  into  (2.58)  results  in 

00      rN     r2ir 

covm(R,T)  =  J2         /      \al{u,0)\2Pl{u,e)e-^Tdwdr,  (2.61) 
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where 


roo     roo 

pi(u,o)  =  J0J0  Mu,*W(<*'**'K{z)C'p(z') 


(2.62) 


rR    rR 


fit 

/      /     Gm(r,  z)Gm{r',  z)  exp[ikhi{to)  cos  0(r  -  r')]drdr' 

Jo    Jo 


dzdz  . 


If  Gm(r,  z)  is  assumed  to  be  range  independent,  then  it  can  be  taken  outside  the 
integral  and  (2.62)  reduces  to 


•R    rR 


pi{u,0)  =   /      /     exp[iku{u>)  cos  0(r  -  r')]drdr  \        ^i{u,z)C'(z)Gm(z)d: 
Jo    Jo  \Jo 


(2.63) 


The  first  integral  on  the  right  hand  side  of  (2.63)  reduces  to 

2[1  —  cos(khi  cos  OR)] 


(khi  cos  oy 


and  can  be  further  reduced  to 


r>2  •      2  ( khiR         n\ 
K  sine    I  — - —  cos  U    , 


where  sinc(x)  =  sinx/x.  Thus  (2.63)  becomes 


-^-cos^J  |yo    Mu,z)C'p(z)Gm(z)d; 


(2.64) 


The  covariance  spectrum  can  be  represented  by 


/oo 
covm{R,  r)e-,QTdr. 
-oo 


(2.65) 


Inserting  (2.61)  into  (2.65)  yields 


00       rN      r2ir  /    roo         .  .  \ 

Fm(R,n)  =  Y,J     Jo     Mu>,0)\  pi(<OiO){J      e^-^TdT)du>de,  (2.66) 


JTiJf 
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which  after  integrating  over  r  gives 


00       rN     flit 


*      fly      rlir 

Fm(R,Q)  =  y£         /     \a,{u>,0)\2p,(u>,0)6(n-u>)du>d0,  (2.67) 

/=i  Jf    Jo 


which  can  be  reduced  to 


00       f2n 


Fm(R,Q)=J2         \a,(n,0)\  pi{n,0)M. 


(2.68) 


The  previous  results  can  be  simplified  by  assuming  that 


\ai(u>t9)\2  =  A(u)B(0)Dh 


(2.69) 


and  further  simplified  by  assuming  B{6)  =  6(6  —  60)  and  that  the  packet  width  of 
internal  waves  is  Rd.  For  the  case  of  6q  =  tt/2,  (2.61)  becomes 


cov, 


fN 
(R)t=o  —  Rd  I 
J  f 


J2Mto)Di 


J%i{il1z)Clp{z)Gm{z)dz 


dil.  (2.70) 


In  this  thesis  we  will  evaluate  equations  (2.53)  and  (2.70).  The  assumptions  we  make 
in  evaluating  those  equations  are  detailed  in  Section  4.3. 
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Chapter  3 

Barents  Sea  Polar  Front 
Experiment 


The  Barents  Sea  Polar  Front  (BSPF)  experiment  took  place  from  August  6  to  August 
26,  1992  in  the  coastal  waters  off  of  Svalbard,  about  42  nautical  miles  due  east  of 
Bj0rn0ya  or  Bear  Island  at  the  northern  edge  of  the  Bj0rn0ya  Trough.  The  BSPF 
experiment  involved  three  CTD  and  ADCP  surveys,  and  the  deployment  of  four 
moorings  for  about  ten  days.  Water  depth  in  the  experimental  area  varied  from 
110  m  to  380  m. 

3.1      Barents  Sea 

The  Barents  Sea  is  a  shallow,  peripheral  basin  bounded  by  Svalbard  on  the  north, 
Norway  and  Russia  on  the  south,  Novaya  Zemlya  on  the  east,  and  the  Norwegian  Sea 
on  the  west.  The  majority  of  the  sea  is  contained  between  70°  and  80°  north  latitude 
and  16°  and  60°  east  longitude  (Fig.  3-1).  The  boundary  between  the  Norwegian 
Sea  and  the  Barents  Sea  is  the  transect  from  Spitsbergen  to  Bear  Island  to  the  coast 
of  Norway.  Table  3.1  gives  the  characteristics  of  the  waters  discussed  below.  For  a 
detailed  description  of  the  water  masses  and  types  see  Hopkins  [1991]  or  Johannessen 
[1986].  The  Barents  Sea  Polar  Front  (BSPF)  is  formed  at  the  boundary  of  two 
currents,  the  North  Cape  Current  to  the  south,  and  the  Bear  Island  Current  to  the 
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Figure  3-1:  Chart  of  the  Barents  Sea  with  major  currents.  The  thick,  solid  line  is  the 
Barents  Sea  Polar  Front. 
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Table  3.1:  Water  Mass  Definitions 
Water  Mass  Symbol         T(°C)  S(ppt)  Comments 

North  Atlantic  Water  NAtW  >  8  >  35.3  Surface  water  of 

N.  Atlantic  Current 
Norwegian  Atlantic  Water     NwAtW  >  2  >  35.0     Water  of  Norwegian 

Atlantic  Current 
Polar  Water  PW     Freezing     33.0-34.4 

to  2 
Barents  Polar  Water  BrPW  -1-3  <  34.3 

north.  The  North  Cape  Current,  an  extension  of  the  Norwegian  Atlantic  Current 
(NwAtC),  carries  warm,  saline  Norwegian  Atlantic  Water  (NwAtW)  into  the  Barents 
Sea.  The  Bear  Island  Current,  which  flows  along  the  northern  part  of  the  Bj0rn0ya 
or  Bear  Island  Trough,  carries  cold,  fresh  Barents  Polar  Water  (BrPW)  out  into  the 
Norwegian  Sea.  The  front  is  a  stable  interface  because  the  warm,  saline  NwAtW 
has  approximately  the  same  density  as  the  cold,  fresh  BrPW.  The  BSPF  is  some 
1500  kilometers  long  and  roughly  follows  the  200  meter  isobath.  It  is  also  strongly 
influenced  by  the  tidal  cycle  and  can  be  displaced  by  a  much  as  10  kilometers  during 
a  tidal  cycle  [Johannessen  1986]. 

The  tides  in  the  Barents  Sea  are  primarily  semidiurnal  (M2)  with  a  weaker  diurnal 
component  (K\).  Both  Gjevik  [1990]  and  Schwiderski  [1986]  have  proposed  tidal 
models  for  the  Barents  Sea.  The  predicted  amplitudes  for  each  model  are  given  in 
table  3.2.  The  models  differ  in  that  Gjevik's  model  resolves  amphidromes  for  both  the 
Mi  and  K\  in  the  region  of  Bear  Island,  whereas  Schwiderski's  does  not.  Huthnance 
[1981]  made  a  series  of  current  meter  measurements  from  Northern  Norway  to  Bear 
Island.  From  his  experiment,  he  estimates  that  both  the  M2  and  K\  are  mostly 
barotropic  (about  86  %  ).  The  predominantly  barotropic  S2  tide  is  smaller  than  the 
M2,  and  has  a  fairly  consistent  relation  with  it. 
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Table  3.2:  Tidal  Amplitudes  at  Bear  Island 


Obs. 

Gjev 

Schw 

Constituent 

Amp. 

Amp. 

Amp. 

(cm) 

(cm) 

(cm) 

M2 

34.2 

34.5 

40 

Kx 

5.4 

6.3 

5.5 

Obs.  —  Observed  values  taken  from  Gjevik  [1990]   . 
Gjev  —  Values  taken  from  Gjevik  [1990] 
Schw  —  Values  taken  from  Schwiderski  [1990] 

3.2  Experimental  Objectives 

The  purpose  of  the  experiment  was  to  study  the  structure  and  the  dynamics  of  the 
Barents  Sea  Polar  Front  by  using  a  variety  of  moored  and  shipboard  oceanographic 
instruments.  A  key  aspect  of  the  BSPF  experiment  was  the  Barents  Sea  Tomography 
Transmission  Test  (BSTTT)  whose  major  goals  were:  first,  to  determine  the  iden- 
tifiability,  resolvability,  strength  and  variability  of  the  multipath  structure  between 
the  acoustic  moorings;  second,  to  examine  the  acoustic  signature  of  various  oceano- 
graphic processes  (surface  waves,  internal  tides,  internal  waves,  eddies,  etc.)  between 
the  acoustic  moorings;  and  third,  to  use  "tomographic  inverse  images"  to  study  the 
region  of  the  BSPF. 

3.3  Data  Acquisition 

Four  acoustic  instruments  were  deployed  during  the  experiment  in  a  trapezoidal  shape 
about  50  kilometers  apart  at  the  base,  with  the  distance  between  the  northeast  and 
southeast  moorings  being  34  kilometers  (Fig.  3-2).  The  Northwest  (NW)  and  South- 
west (SW)  moorings  contained  400  Hz  tomography  transceivers  while  the  Northeast 
(NE)  mooring  contained  a  224  Hz  projector.  These  three  moorings  were  deployed  for 
a  period  of  9  days  (NE)  to  12  days  (SW).  The  transceiver  on  the  NW  mooring  acted 
as  400Hz  receiver  only,  while  the  transceiver  on  the  SW  mooring  acted  intermittently 
as  a  400Hz  source  and  continuously  as  a  400Hz  receiver.    The  fourth  or  Southeast 
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Figure  3-2:  Chart  of  the  BSPF  experiment.  Triangles  are  array  positions.  Arrows 
denote  the  acoustic  paths  during  the  experiment. 

(SE)  mooring  had  a  vertical  line  array(VLA)  of  sixteen  hydrophones  which  collected 
70  hours  of  data.  The  data  from  the  VLA  was  telemetered  in  real  time  back  to 
the  research  ship  using  a  wireless  local  area  network.  The  acoustic  paths  between 
the  moorings  are  shown  in  Fig.  3-2.  Details  of  the  vertical  line  array  are  given  in 
von  der  Heydt  et  al.  [1992].  All  of  the  moorings  except  the  SW  mooring  contained 
several  current  meters  and  Brancker  temperature  loggers  to  provide  a  vertical  profile 
of  current  and  temperature  at  each  mooring.  Details  are  given  in  Table  3.3. 

Mooring  motion  of  the  source  or  receiver  due  to  tides  or  other  currents  changes 
the  source-receiver  range,  and  introduces  error  into  the  measurement  of  the  travel 
time.  To  correct  for  this  mooring  motion,  all  of  the  moorings  except  for  the  NE 
mooring  used  acoustic  navigation  systems.  The  determination  of  mooring  motion  is 
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Figure  3-3:  CTD  grid  employed  during  BSPF  experiment.  Circles  are  CTD  stations. 
The  grid  was  fully  covered  three  times. 

accurate  to  within  lm.  Additionally,  the  SE  mooring  had  two  tilt  meters.  The  ship 
also  conducted  three  complete  70  km  by  80  km  CTD  grids  (Fig.  3-3)  of  the  research 
area,  two  detailed  grids  across  the  Barents  sea  polar  front,  and  a  24  hr  CTD  anchor 
station  at  the  SW  mooring.  Details  of  the  CTD  grids  and  the  anchor  station  are  given 
in  Pawlowicz  [1992].  Several  SUS  charges  were  also  dropped  to  allow  determination 
of  bottom  properties  and  reverberation  characteristics. 
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Table  3.3:  Barents  Sea  Moorings 

Instrument  Depth  (m) 


Northwest  Mooring— 170m 


Northeast  Mooring-142m 


Southwest  Mooring-380m 


Southeast  Mooring-280m 


Neil  Brown  SACM 
Brancker  Tlog 
Brancker  Tlog 
Aanderaa  CM 
400  Hz  transc. 
Aanderaa  CM 


Neil  Brown  ACM 
Brancker  Tlog 
Brancker  Tlog 
Neil  Brown  ACM 
Aanderaa  CM 
224  Hz  proj. 


Neil  Brown  SACM 
Brancker  Tlog 
Brancker  Tlog 
Aanderaa  CM 
Aanderaa  CM 
Brancker  Tlog 
Brancker  Tlog 
400  Hz  transc. 


Vertical  Array 


19.4 
29.5 
39.2 
49.3 
66.7 
159.9 


22.1 

32.1 
42.1 
52.0 
82.1 
121.8 


19.9 

29.9 

39.9 

50.0 

80.0 

140.0 

200.5 

358.1 


126-275.8 
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Chapter  4 


Temperature  Field  Data  Analysis 


The  raw  data  from  which  the  calculations  in  this  chapter  are  made  comes  mainly  from 
the  CTD  casts  made  at  each  mooring  site,  and  the  moored  Brancker  temperature 
loggers.  Some  temperature  information  is  also  from  sensors  in  the  current  meters. 
Figure  4-1  is  a  typical  Brancker  temperature  logger  record.  The  temperature  record 
is  non-linear  in  that  at  higher  temperatures  the  signal  is  "spikier"  than  at  lower 
temperatures.  This  effect  is  most  pronounced  for  those  loggers  located  where  the 
temperature  gradient  is  most  rapidly  changing.  A  one  meter  downward  displacement 
of  water  at  the  "knee"  of  the  temperature  curve  (i.e.  the  mixed  layer  base)  results  in 
a  much  greater  change  of  temperature  than  a  one  meter  upward  displacement.  Figure 
4-2(a)  shows  a  sinusoidal  9.9m  amplitude  particle  displacement  which,  acting  on  a 
typical  temperature  profile  (Fig.  4-2(b))  results  in  a  "distorted"  temperature  record 
(Fig.  4.2(c)). 

The  changes  in  temperature  seen  are  due  to  both  a  vertical  particle  displacement 
of  water  caused  by  internal  tides  and  waves,  and  a  horizontal  advection  of  water 
caused  by  the  barotropic  tide.  Which  mechanism  is  dominant  can  be  determined 
by  evaluating  constraints  imposed  by  the  local  temperature  gradient,  current  veloc- 
ities, and  water  depth.  Figure  4-3  shows  three  CTD  temperature  profiles  taken  at 
the  Southwest  array.  If  we  take  profile  #418  as  the  baseline,  and  translate  profile 
#218  vertically  by  +7m,  then  profile  #218  will  lie  directly  over  profile  #418.  This 
corresponds  to  a  temperature  change  solely  due  to  internal  tides  and  internal  waves. 
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Figure  4-1:  Temperature  record  from  Brancker  logger  #  3666  located  on  the  SW 
mooring  at  a  depth  of  29.9m.  Brancker  temperature  resolution  is  0.1  °C 

Similarly,  if  profile  #118  is  translated  upwards  in  depth  by  15m  and  moved  hori- 
zontally by  0.25  °C  it  will  lie  over  profile  #418.  This  corresponds  to  a  temperature 
change  due  to  both  a  barotropic  current  and  internal  tides  and  internal  waves.  We  will 
now  consider  the  different  effects  separately,  and  thereby  determine  the  percentage 
that  each  signal  contributes  to  the  whole. 

4.1      Horizontal  Advection  versus  Internal  Tides 
and  Waves 

The  amount  of  horizontal  advection  was  determined  at  each  of  the  three  instrumented 
arrays  by  using  the  assumption  that  the  temperature  changes  were  solely  due  to 
horizontal  advection  and  that  any  vertical  motion  due  to  internal  tides  and  waves 
or  topographic  changes  was  negligible.  In  order  to  explain  the  method,  the  data 
from  the  NE  mooring  will  be  used.    Figure  4-4  shows  the  temperature  and  velocity 
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Figure  4-2:  (a)  is  a  9.9m  sinusoidal  particle  displacement  about  the  logger  depth 
(dashed  line),  (b)  is  a  smoothed  temperature  profile  at  the  SW  mooring.  The  dashed 
line  is  the  logger  depth,  (c)  is  the  "distorted"  temperature  record  which  results  from 
the  sinusoidal  particle  displacement  in  a  varying  temperature  gradient.  The  dashed 
line  is  the  axis  of  zero  phase. 
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Figure  4-3:  Three  temperature  profiles  at  the  SW  array  taken  over  a  span  of  18 
days.  The  temperature  records  show  both  a  temperature  change  due  to  vertical 
displacement  (caused  by  internal  tides  and  waves)  and  due  to  horizontal  advection 
(caused  by  the  barotropic  tide) 

records  from  Aanderaa  #9259  at  a  depth  of  82.1  meters.  The  temperature  record  has 
a  minimum  of  —0.9  °C  and  a  maximum  of  —0.36  °C  with  a  peak  to  peak  change  in 
temperature  of  0.208  °C. 

Figure  4-5  is  a  contour  map  of  temperatures  taken  from  the  CTD  profiles  sur- 
rounding the  NE  mooring  and  shows  that  the  warmer  water  at  82  meters  is  to  the 
east  and  to  the  south,  with  the  east-west  gradient  being  much  larger  than  the  north- 
south  gradient.  Several  things  are  apparent  in  Fig.  4-4.  The  u  velocity  lags  the 
temperature  record  by  90  degrees,  while  the  v  velocity  is  in  phase  with  the  tempera- 
ture record.  Given  that  the  u  velocity  is  over  twice  the  v  velocity  and  given  that  the 
east  to  west  temperature  gradient  is  much  larger  than  the  north  south  temperature 
gradient,  the  predominant  effect  will  be  from  the  u  velocity.  Since  the  u  velocity  lags 
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Figure  4-4:  U  and  V  velocity,  and  Temperature  records  from  Aanderaa  #9259  at  the 
NE  array  at  a  depth  of  82.1  meters. 
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Figure  4-5:  Temperature  contour  map  at  NE  mooring  at  a  depth  of  82.1m  using 
adjacent  CTD  stations.  The  current  ellipse  from  Aanderaa  #  9259  shows  the  amount 
of  horizontal  advection  during  one  tidal  cycle. 
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temperature  by  90  degrees,  at  the  time  when  the  current  shifts  from  flowing  westward 
to  flowing  eastward,  the  temperature  will  be  at  a  maximum,  i.e.  it  is  at  this  time  that 
the  warm  water  to  the  east  will  have  been  advected  the  greatest  distance  to  the  west. 
This  can  also  be  seen  by  taking  the  integral  of  the  u  velocity;  the  displacement  record 
is  now  180  degrees  out  of  phase  with  the  temperature  record  which  is  consistent  with 
having  the  warmer  water  to  the  east. 

The  displacement  ellipse  was  then  plotted  over  the  temperature  contour.  From  this 
plot,  one  sees  that  the  maximum  change  in  temperature  due  to  horizontal  advection 
is  0.037  °C,  which  only  accounts  for  18  percent  of  the  peak  to  peak  temperature 
change.  Another  argument  against  horizontal  advection  being  the  mechanism  behind 
the  temperature  changes  is  that  the  —0.9  °C  water  would  have  had  to  been  advected 
from  over  14  km  away  from  the  northwest  and  the  —0.36  °C  water  would  have  had  to 
been  advected  from  over  10  km  away  from  the  east  or  over  20  km  away  from  the  south. 
None  of  these  distances  are  credible  given  the  results  listed  in  Table  4.1.  Analyses  of 

Table  4.1:  Horizontal  Advection  of  Water  due  to  Barotropic  Tides 

Array  Depth  (m)  N-S  disp  (km)  E-W  disp  (km) 

NE              82J                   6\98  2L2 

NW             49.3                   1.76  2.46 

SW              80.0                   0.89  1.82 

the  S W  and  N W  moorings  yield  similar  results. 

4.2      The  Internal  Tide  and  Internal  Wave  Spec- 
trum 

In  this  section  we  convert  temperature  fluctuations  into  particle  displacement  fluctua- 
tions and  then  into  internal  tide  and  internal  wave  amplitudes  using  Fourier  analysis. 
In  performing  this  analysis  we  assumed  that  the  temperature  fluctuations  can  be 
used  to  obtain  the  magnitude  of  the  particle  displacement.  This  will  be  true  if  the 
temperature-salinity  (T-S)  relationship  is  constant,  i.e.  one  can  go  uniquely  from 
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Figure  4-6:  Temperature  and  Salinity  versus  depth  for  the  CTD  casts  taken  at  the 
NE  array. 

temperature  displacement  to  density  displacement.  Figure  4-6  shows  the  relative 
constancy  of  the  T-S  relationship  over  three  different  CTD  casts  at  one  site.  A  rise 
in  temperature  for  a  sensor  (see  Fig.  4-7)  above  the  mixed  layer  (i.e.  negative  tem- 
perature gradient)  corresponds  to  a  downward  displacement  of  warmer  water.  The 
vertical  displacement  6z  for  a  parcel  of  water  is 


dz 

6z  =  6T—, 
dT 


(4.1) 


where  ST  is  the  temperature  fluctuation  measured  by  the  Brancker  logger,  and  dz/dT 
is  the  temperature  gradient  at  the  mooring.  We  have  taken  z  as  being  negative 
downward.  The  particle  displacement  can  be  determined  by  integrating  (4.1)  so  that 
each  temperature  corresponds  to  a  given  depth.  This  method  fails  if  the  temperature 
profile  is  double  valued  near  the  depth  of  the  logger  or  the  logger  is  near  the  knee 
of  the  T(z)  profile.  Using  the  temperature  profiles  from  the  CTD  casts,  we  obtained 
a  best  fit  temperature  profile  which  preserved  the  slopes  of  the  CTD  temperature 
profiles.     Each  temperature  record  was  then  converted  to  a  particle  displacement 
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Figure  4-7:  Temperature  record  for  Brancker  logger  #3663  at  NE  mooring  at  a  depth 
of  32.1m. 

record  (see  Fig.  4-8).  The  method  assumed  that  the  best-fit  linearized  temperature 
profile  was  temporally  stationary,  whereas  in  fact  it  is  not  quite,  as  demonstrated 
by  Fig.  4-3.  Additionally,  in  several  records  the  temperature  minima  and  maxima 
recorded  by  the  Brancker  temperature  loggers  exceeded  the  temperature  minima  and 
maxima  recorded  during  the  CTD  casts.  The  CTD  casts  represent  "snapshots"  of  the 
temperature  profile,  and  of  the  ten  casts  taken  at  the  NE,  SW,  or  NW  moorings  only 
two  occurred  (one  at  the  NW  mooring  and  one  at  the  NE  mooring)  when  the  moorings 
were  deployed,  and  these  CTD  casts  were  at  the  ends  of  the  temperature  records.  In 
order  to  compensate  for  this  mismatch,  the  temperature  profiles  were  extended  at  the 
upper  and  lower  ends  while  maintaining  the  same  slope.  This  process  was  analogous 
to  including  some  horizontal  displacement.  It  must  be  remembered  that  there  are 
two  processes  taking  place,  horizontal  advection  and  vertical  displacement,  and  that 
when  the  Brancker  temperature  exceeds  the  CTD  temperature,  what  has  occurred  is 
obviously  both  a  horizontal  advection  and  a  vertical  temperature  displacement. 

To  estimate  the  error  of  extending  the  temperature  profiles,   the  temperature 
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Figure  4-8:    Particle  displacement  (  calculated  for  Brancker  logger  #  3663  on  NE 
array  at  a  depth  of  32.1m. 


record  for  Brancker  #3660  was  examined.  This  temperature  record,  like  the  others  at 
the  SW  mooring,  showed  several  eddies,  which  act  similar  to  a  horizontal  advection  of 
the  profile  (e.g.  Fig.  (4-1).  The  particle  displacement  and  amplitudes  were  computed 
two  ways  for  the  Brancker  #3660  temperature  record:  first,  assuming  that  the  eddy 
was  a  horizontal  advection  and  adjusting  the  temperature  profile  accordingly;  and 
second,  filtering  out  the  eddy  feature  and  using  an  unadjusted  temperature  profile. 
The  difference  between  the  amplitudes  in  the  internal  tide  band  computed  the  two 
different  ways  was  at  most  less  than  3.5%. 

Once  the  particle  displacement  was  computed,  it  was  then  Hanning  windowed,  and 
analyzed  using  the  Welch  method  of  fast  Fourier  transform(FFT)  spectrum  analysis. 
Figure  4-9  is  a  typical  spectrum.  The  most  significant  peak  in  the  internal  tide  region 
is  the  semidiurnal  peak.  The  amplitude  was  extracted  by  summing  the  area  under 
the  peak 


A  = 


2  k=b 

jffpm 


(4.2) 
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Figure  4-9:  Power  spectrum  density  for  the  particle  displacement  calculated  for 
Brancker  logger  #  3663  on  the  NE  array  at  a  depth  of  32.1  m.  The  estimate  has 
four  degrees  of  freedom. 

where  A  is  the  amplitude,  PD(k)  is  the  power  density  function  at  k  and  TV  is  the  length 
of  the  sequence.  Table  4.2  gives  the  amplitudes  at  the  Southwest  array  computed 
using  (4.2)  for  the  M2  tidal  component. 

Due  to  the  shortness  of  the  records,  high  resolution  Fourier  analysis  is  not  possible. 
This  results  in  combined  tidal  amplitudes  for  some  frequency  values.  The  M2  peak 
probably  contains  52  and  other  semidiurnal  constituents,  though  they  are  small.  Also, 
the  amplitudes  computed  this  way  are  average  amplitudes,  and  therefore  smooth  out 

Table  4.2:  SW  Mooring  Particle  Displacement  Amplitudes  £ 
Logger  Depth(m)     M2  (m) 
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Figure  4-10:  Amplitude  spectrum  for  Brancker  logger  #  3660  on  SW  array  at  a  depth 
of  39.9m. 

the  previously  noted  intermittent  nature  of  the  internal  tides  and  internal  waves. 
This  method  gives  reliable  estimates  for  frequencies  in  the  internal  tide  range  (1-10 
cpd),  but  may  underestimate  the  amplitudes  in  the  higher  frequency  internal  wave 
band.  This  effect  will  be  addressed  when  estimating  the  amplitude  spectrum  used  in 
equation  (2.70).  Figure  4-10  shows  a  typical  amplitude  spectrum. 

The  amplitudes  extracted  using  Fourier  analysis  are  used  to  determine  the  relative 
weights  of  each  of  the  normalized  internal  tide  and  internal  wave  modes.  This  analysis 
is  performed  for  the  SW  mooring  because  it  contains  four  Brancker  temperature 
loggers  and  spans  170m  of  the  water  column,  while  the  NE  and  NW  moorings  only 
have  two  temperature  loggers  and  span  only  10m  of  the  water  column.  Also,  using 
a  least  squares  fit,  up  to  four  modes  can  be  estimated  at  the  Southwest  array  (Fig. 
4- 11  (a)).  The  least  squares  fit  was  simplified  by  averaging  together  the  two  shallow 
loggers,  as  they  are  only  10m  apart  and  can  be  approximated  as  a  single  sample. 
Figure  4- 11(b)  is  the  result  of  this  average.  The  results  from  both  methods  were 
consistent.    The  fit  was  only  performed  for  frequencies  up  to  30  cpd,  since  this  is 
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the  maximum  buoyancy  frequency  for  the  two  deep  loggers,  and  in  fact,  the  method 
appears  to  break  down  at  about  27cpd.  Also,  the  sign  of  the  mode  shapes  was 
consistent  in  the  internal  tide  band  (0-8  cpd),  but  the  sign  of  the  third  mode  changed 
in  the  frequency  band  of  10-30  cpd.  This  change  coincided  with  change  from  the 
logger  at  140m  having  the  largest  amplitude  in  the  internal  tide  bank  to  the  logger  at 
200m  having  the  largest  amplitude  in  the  internal  wave  band.  The  percentages  were 
averaged  over  the  internal  tide  and  internal  waves,  resulting  in  an  average  percentage 
for  each  of  the  first  three  modes  of:  mode  #1  -  0.58,  mode  #2  -  0.26,  and  mode  #3 
-  0.16.  As  a  check  on  these  numbers,  the  modes  were  also  estimated  incoherently,  by 
performing  a  least  squares  fit  of  the  first  mode  squared  to  the  square  of  the  amplitudes. 
The  residual  was  then  least  square  fitted  to  the  second  mode  squared.  This  process 
was  performed  for  the  first  three  modes  at  the  Ml  tidal  frequency.  The  resulting 
relative  percentages  for  the  first  three  modes  were:  mode  #1  -  0.62,  mode  #2  -  0.22, 
and  mode  #3  -  0.16.  Both  of  these  sets  of  numbers  are  consistent  with  each  other 
and  the  results  reviewed  in  Wunsch  [1975]. 

4.3      Review  of  Assumptions 

We  are  now  ready  to  compute  the  acoustic  fluctuations  using  Eqs.  (2.53)  and  (2.70), 
but  before  doing  so  a  review  of  the  validity  of  the  assumptions  is  warranted.  First, 
we  have  assumed,  in  truncating  the  IW  mode  series  to  three  modes,  that  the  internal 
tide  and  internal  wave  spectra  are  red  i.e.  most  of  the  energy  is  in  the  lowest  modes. 
This  is  probably  an  excellent  assumption,  as  numerous  experiments  have  concluded 
that  the  bulk  of  the  energy  in  both  the  internal  tide  and  internal  wave  bands  is  in  the 
lowest  modes.  The  ratios  of  the  first  three  internal  wave  modes  calculated  in  Section 
4.2  show  good  agreement  with  those  in  Wunsch  [1975]. 

Second,  we  have  assumed  in  solving  equation  (2.70)  that  the  direction  of  propaga- 
tion is  along-shelf  instead  of  across-shelf.  We  assumed  this  direction  of  propagation 
for  two  reasons:  first,  we  didn't  have  available  with  our  single  point  moorings  any 
information  on  the  direction  of  propagation;  and  second,  this  direction  (90  =  7r/2) 
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Figure  4-11:  (a)  is  the  least  squares  fit  of  four  amplitudes  measured  at  the  SW  array 
to  a  sum  of  the  first  four  IW  modes,  (b)  is  the  same  calculation  except  that  the 
two  shallow  Branckers  are  averaged  together.  In  both  figures,  considerable  noise  is 
evident  in  the  internal  tide  region,  and  the  method  is  unreliable  above  about  27cpd 
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maximizes  the  sine2  function  and  thereby,  the  travel  time  perturbation  covariance 
(Eq.  2.61).  This  gives  the  largest  effect  we  could  hope  to  see,  and  thus  represents 
an  "upper  limit"  calculation.  Additionally,  although  the  assumption  of  along-shelf 
propagation  is  more  tenuous  than  that  of  across-shelf  propagation,  it  is  perhaps  not 
wholly  unjustified,  since  the  experimental  region  is  on  the  edge  of  the  Bear  Island 
trough  and  in  the  vicinity  of  the  strong  Barents  Sea  Polar  front.  As  discussed  in 
Section  2.2.2,  although  internal  tides  and  internal  waves  propagate  shoreward,  one 
experiment  noted  that  the  internal  tides  propagated  along  the  coast.  It  was  postu- 
lated that  either  coastal  fronts,  or  mesoscale  structures  were  refracting  the  internal 
tide  or  that  internal  tides  radiating  from  differently  oriented  shelves  were  causing  a 
complicated  wave  structure  [Levine  and  Richman,  1989].  The  experimental  area  con- 
tains both  such  features.  The  Bear  Island  trough  certainly  qualifies  as  a  shelf  area 
with  different  orientations. 

Third,  we  have  assumed  that  the  change  in  temperature  seen  by  the  logger  was 
largely  due  to  a  vertical  particle  displacement  caused  by  internal  tides  and  internal 
waves  rather  than  due  to  a  horizontal  advection  caused  by  the  barotropic  tide.  We 
showed  in  Section  4.1  that  at  most  18%  of  the  change  in  temperature  could  be  at- 
tributed to  horizontal  advection.  This  effect  was  included  implicitly  in  determining 
the  particle  displacements,  but  by  not  including  it  resulted  in  a  less  than  3.5%  error. 

Fourth,  we  have  assumed  that  <SC(x,  t)  ^>  u(x,  t)  •  r.  From  the  Aanderaa  data 
the  maximum  current  was  0.325  m/s,  whereas  in  the  mixed  layer,  a  5m  particle 
displacement  could  result  in  a  6.0  m/s  change  in  sound  speed. 
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Chapter  5 

Travel  Time  Disturbance 
Predictions 


Using  the  theory  from  chapter  two,  and  the  internal  wave  frequency  spectrum  and 
amplitudes  generated  using  the  methods  of  chapter  four,  a  modal  picture  and  a  ray 
picture  of  the  travel  time  perturbations  are  now  generated.  These  predictions  will 
eventually  be  compared  with  the  actual  results  from  the  experiment,  when  these 
become  available  from  colleagues  of  the  Naval  Postgraduate  School  who  are  working 
on  the  travel  time  estimation  problem. 

5.1      Theoretical  Results  for  the  M2  Tide 

In  chapter  two,  a  formulation  for  the  travel  time  perturbation,  equation  (2.53),  and 
travel  time  covariance,  equation  (2.70),  was  derived,  using  adiabatic  normal  modes  as 
a  basis.  The  adiabatic  assumption  requires  a  "slowly"  varying  sound  speed  profile  and 
bathymetry.  We  further  assumed  range  independence  for  the  group  velocity  dispersion 
Gm(z)  and  for  the  particle  displacement  (.  In  evaluating  equations  (2.53)  and  (2.70), 
we  chose  a  simple,  two-segment,  acoustic  and  oceanographic  model.  The  model  was 
generated  by  dividing  the  acoustic  path  from  the  Northeast  mooring  to  the  Southeast 
mooring  into  two,  flat-bottomed  segments.  Segment  number  1  is  north  of  the  BSPF, 
while  segment  number  2  is  south  of  the  front.   Each  segment's  depth  is  the  average 
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Figure  5-1:  Bathymetry  of  two  segment  range  dependent  model  from  NE  array  to  SE 
array.  Solid  line  is  the  depth  of  the  model.  Dashed  line  is  the  depth  at  each  CTD 
station.  Segment  #1  goes  from  0-26. 93km. 

of  the  depths  recorded  at  each  of  the  CTD  stations  between  the  moorings.  Figure 
5.1  shows  the  bathymetry  (using  CTD  station  depths)  superimposed  on  the  segment 
depths.  The  physical  properties  (temperature,  salinity,  et  al.)  are  the  averages  from 
CTD  stations  between  the  moorings.  Figures  5-2(a)  and  (b)  show  the  average  sound 
profiles  of  each  segment. 

The  total  travel  time  disturbance  will  be  taken  as  the  sum  of  the  disturbances  in 
the  two  segments.  This  is  somewhat  inconsistent  in  that  it  ignores  mode  coupling 
effects  due  to  the  front.  However,  the  exact  magnitude  of  the  mode  coupling  due  to 
the  front  is  not  well  known  at  this  time,  so  we  did  not  pursue  this  correction  to  our 
result. 

5.1.1      Acoustic  Ray  Results 

The  simplest,  shallow  ocean  model,  given  a  strong  mixed  layer,  is  a  two-layer  ocean. 
In  this  model,  the  sound  speed  of  each  layer  is  the  vertical  average  of  the  Fig.   5-2 
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Figure  5-2:    (a)  Averaged  sound  speed  profile  for  segment  #  1,  and  (b)  Averaged 
sound  speed  profile  for  segment  #  2 

sound  speed  profile  for  a  given  range  segment.    Table  5.1  lists  the  details  for  each 
segment. 

The  travel  time  perturbations  in  ray  theory  can  easily  be  calculated  by  perturbing 
the  layer  an  amount  equal  to  the  peak  to  peak  particle  displacement  £  and  inputting 
the  resultant  sound  speed  perturbation  8C  into  a  discretized  form  of  equation  (2.7). 
This  gives 

(5.1) 


ST=T  26C<<z^NAz 


where  N  is  the  number  of  complete  ray  cycles,  and  \  is  the  grazing  angle. 

This  calculation  was  run  on  segment  #  1  for  the  M2  internal  tide  using  a  peak- 
to-peak  particle  displacement  £  of  11  m  with  a  source  and  receiver  in  the  lower  layer. 
It  was  assumed  to  first  order  that  neither  \  nor  N  changed,  that  is  the  ray  path  was 
identical  before  and  after  the  sound  speed  profile  was  perturbed.  Figure  5-3(a)  shows 
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the  results  for  eigenrays  which  penetrate  the  upper  (mixed)  layer. 

Figure  5-3(b)  is  the  calculation  performed  for  segment  #  2  with  a  peak-to-peak 
particle  displacement  £  of  15  m,  which  was  the  maximum  observed  displacement  of 
the  thermocline  seen  in  Fig.  4-3.  The  displacement  of  the  thermocline  due  to  the 
M2  internal  tide  may  be  as  small  as  4  m,  which  was  the  computed  Ml  internal  tide 
amplitude  at  the  SW  array  using  Fourier  analysis.  Using  the  smaller  value  of  4  m 
would  scale  Fig.  5-3(b)  by  one  fourth. 

The  calculation  can  also  be  performed  using  a  more  sophisticated  ray  tracing 
program  such  as  MPP  with  fully  range  dependent  profiles  and  geometry.  As  we  are 
only  interested  in  order  of  magnitude  estimates  with  this  ray  calculation,  we  have  not 
pursued  this. 

5.1.2      Acoustic  Normal  Mode  Results 

In  calculating  the  acoustic  modal  results,  the  same  two-segment  view  of  the  acoustic 
path  was  taken.  Figures  5-4(a)  and  (b)  show  buoyancy  frequency  profiles  calculated 
from  equation  (2.25)  using  averaged  temperature  and  salinity  profiles  for  segments  # 
1  and  #2. 

Using  these  buoyancy  frequency  profiles,  equation  (2.27)  can  be  solved  to  give  the 
internal  tide  and  internal  wave  modes.  Figures  5-5(a)  and  (b)  show  the  first  three 
internal  tide  modes  for  each  segment  at  the  M2  frequency.  Similarly,  given  the  sound 
speed  profiles,  equation  (2.41)  gives  the  adiabatic  normal  modes  (Figs.  5-6(a)  and 
(b)).  In  comparing  Fig.  5-5  with  Fig.  5-6,  the  case  can  be  made  that  the  third  Ml 
internal  tide  mode  will  possibly  have  the  greatest  effect  on  the  acoustics,  because  its 

Table  5.1:  Two-layer,  range  independent  model 

Segment  Nr.  1  Segment  Nr.  2 

Distance  (m)                                           26,930  7680 

Avg.  water  depth  (m)                              185  260 

Mixed  layer  depth  (m)                               25  25 

Sound  speed  (m/s)  upper  layer              1473  1479.9 

Sound  speed  (m/s)  lower  layer               1453  1465.1 
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Figure  5-3:  Travel  time  perturbation  of  eigenrays  versus  launch  angle. 
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Figure  5-4:    Smoothed  buoyancy  frequency  N  profiles  for  (a)  segment  #  1  and  (b) 
segment  #  2 
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Figure  5-5:  First  three  internal  tide  modes  for  (a)  segment  #  1  and  (b)  segment  #  2 
at  u  =  Ml 
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Figure  5-6:  First  three  adiabatic  normal  modes  for  (a)  segment  #  1  and  (b)  segment 

#2. 

maximum  occurs  near  the  maximum  of  the  acoustic  modes. 

The  results  of  calculations  using  equation  (2.53)  are  shown  in  Figs.  5-7  to  5-9 
in  three  different  ways.  First,  as  shown  in  Fig.  5-7,  we  assume  that  the  lowest  ten 
internal  tide  modes  have  an  equal  amplitude  of  5.2  m.  This  calculation  allows  us 
to  see  how  the  internal  tide  and  acoustic  modes  will  interact,  everything  else  being 
equal.  For  segment  #1,  the  maximum  travel  time  perturbation,  indeed,  occurs  for 
the  third  internal  tide  mode.  Except  for  the  first  internal  tide  mode,  low  internal 
tide  modes  dominate  the  scattering.  The  internal  tide  modes  have  little  effect  on 
the  first  acoustic  mode.  For  segment  #2,  the  low  internal  tide  modes,  including  the 
first  internal  tide  mode,  also  dominate  the  scattering,  but  the  first  acoustic  mode  is 
somewhat  more  affected  by  the  internal  tides.  It  should  be  noted  that  one  reason 
that  the  travel  time  perturbations  for  segment  #  2  are  much  smaller  than  those  for 
segment  #  1  is  that  the  range  for  the  former  is  four  times  smaller  than  the  latter. 

For  the  second  method,  shown  in  Fig.  5-8,  only  the  three  lowest  internal  tide 
modes  are  used,  instead  of  using  the  lowest  ten  modes  as  before.  The  three  internal 
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Figure  5-7:  Travel  time  perturbations  at  A/2  using  first  ten  internal  tide  modes 
equally  weighted  with  a  maximum  amplitude  of  5.2m.  Lines  are  acoustic  modes: 
#l(o),  #2(+),  #3(x),  #4(*),  #5(.). 
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tide  modes  are  no  longer  equally  weighted,  but  are  now  weighted  by  the  amplitude 
factors,  0.58,  0.26,  and  0.16,  respectively,  which  were  calculated  in  Section  4.2.  This 
calculation  allows  us  to  see  the  contribution  of  each  internal  tide  mode  to  the  acoustic 
travel  time  perturbation.  In  segment  #1,  we  see  again  that  first  internal  tide  mode 
has  little  effect  on  the  low  acoustic  modes,  and  that  the  low  internal  tide  modes  (mode 
two  and  three)  dominate  the  low  acoustic  modes.  In  segment  #2  the  first  internal 
tide  mode  is  dominant  for  all  the  acoustic  modes  calculated. 

For  the  third  method,  the  three  lowest  internal  tide  modes  are  coherently  summed 
together  using  the  signs  calculated  in  Section  4.2.  The  result  is  the  solid  line  in  Fig. 
5-9.  The  dashed  line  in  Fig.  5-9  is  the  contribution  from  the  first  internal  tide  mode 
(the  solid  line  in  Fig.  5-8).  In  both  segments  #1  and  #2,  the  first  mode  is  the 
dominant  contributor  to  the  travel  time  perturbation.  For  segment  #1  (Fig.  5-8)  the 
magnitudes  of  the  second  and  third  internal  tide  modes  are  greater  than  that  of  the 
first  internal  tide  mode  over  the  range  of  the  first  four  acoustic  modes,  but  because 
the  signs  of  the  second  and  third  internal  tide  modes  are  opposite  at  depths  below  60 
m  (where  the  group  velocity  dispersion  is  a  maximum  at  the  lower  acoustic  modes), 
the  effects  of  the  second  and  third  internal  tide  modes  largely  cancel  each  other  out, 
thus  making  the  first  internal  tide  mode  the  largest  component  of  the  travel  time 
perturbation.  For  segment  #2  this  effect  is  less  pronounced,  because  the  magnitudes 
of  the  second  and  third  internal  tide  modes  are  dissimilar.  Figure  5-10  gives  the 
root  mean  square  of  the  travel  time  perturbations,  the  square  root  of  the  solution  to 
equation  (2.70),  for  segments  #1  and  #2  using  the  first  three  internal  tide  modes. 
The  rms  of  the  variance  of  the  travel  time  perturbations  for  both  segments  is  only 
slightly  smaller  than  the  respective  travel  time  perturbations. 

5.1.3      Comparison  of  Ray  and  Modal  Results 

For  segment  #  1  the  Fig.  5-9(a)  modal  travel  time  perturbations  of  0.5  to  35.5  msecs 
are  of  the  same  order  as  the  ray  travel  time  perturbations  of  7.6  to  21.3  msecs.  For 
segment  #  2  the  Fig.  5-9(b)  modal  travel  time  perturbations  of  0.7  to  6.4  msecs  are  of 
the  same  order  as  the  ray  travel  time  perturbations  of  2.7  to  4.1  msecs.  This  gives  us 
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Figure  5-8:    Travel  time  perturbations  at  M2.    The  first  3  internal  tide  modes  are 
shown  and  are  weighted  by  0.58,  0.26,  and  0.16,  respectively. 
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Figure  5-9:  Travel  time  perturbations  at  M2  using  the  sum  of  the  first  3  internal  tide 
modes,  weighted  by  0.58,  0.26,  and  0.16.  Shown  for  comparison  (as  a  dashed  line)  is 
the  contribution  from  the  first  internal  tide  mode  at  M2  (Fig.  5-8). 
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Figure  5-10:  RMS  of  the  travel  time  perturbations  at  M2  using  the  first  three  internal 
tide  modes. 
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some  confidence  that  our  calculations  are  at  least  consistent,  within  the  assumptions 
we  have  made. 

5.2      Acoustic  Normal  Mode  Results  for  the  Inter- 
nal Tide  and  Internal  Wave  Fields 

Section  5.1.2  evaluated  the  travel  time  perturbations  and  covariances  due  to  the 
largest  internal  tide  component,  the  M2  tide.  In  this  section,  using  the  same  two 
segment  model  we  perform  the  same  calculations  over  the  internal  tide  and  internal 
wave  band.  Figures  5- 11  (a)  and  (b)  are  the  travel  time  perturbations  over  the  entire 
IW  band  for  segments  #  1  and  #2.  Although  the  travel  time  perturbation  over  the 
entire  internal  wave  band  has  a  similar  shape  to  the  travel  time  perturbation  from 
the  first  IW  mode  at  the  M2  tidal  frequency  (dashed  line  in  Fig.  5-11),  the  internal 
wave  band  contributes  significantly  to  the  travel  time  perturbation.  Modal  travel 
time  covariances  were  then  computed  using  equation  (2.70).  The  lowest  ten  acoustic 
modes,  and  the  lowest  3  internal  tide  and  internal  wave  modes  were  used  in  the 
calculations.  The  equation  was  discretized  and  summed  from  the  inertial  frequency 
/  to  the  buoyancy  frequency  TV.  The  calculation  is  a  conservative  estimate,  in  that 
in  the  internal  tide  band,  the  tidal  components  were  treated  as  delta  functions,  and 
only  the  A/2  tidal  component  was  included.  The  4  cpd  and  8  cpd  components  were 
omitted,  as  their  significance  is  doubtful  (see  Fig.  4-9).  The  internal  wave  band 
from  10  cpd  to  the  buoyancy  frequency  N  was  treated  as  a  continuous  spectrum,  and 
numerically  integrated.  Figure  5-12  gives  the  results  for  the  acoustic  modes.  The  rms 
of  the  variance  of  the  travel  time  perturbations  for  both  segments  has  only  increased 
slightly,  when  the  higher  frequencies  in  the  internal  wave  band  are  included.  This  is 
because  the  internal  wave  band  is  an  order  of  magnitude  down  compared  to  the  M2 
tide  in  spectral  amplitude. 

When  we  look  at  directional  effects  on  the  internal  waves  and  tides,  the  dominance 
of  the  internal  tide  signal  over  the  internal  wave  signal  becomes  even  greater.  At  the 
M2  tidal  frequency,  scattering  due  to  waves  in  the  across  shelf  direction  (0  =  0)  is  only 
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Figure  5-11:  Travel  time  perturbations  using  the  first  three  IW  modes  from  /  to  N. 
They  are  weighted  by  0.58,  0.26,  and  0.16.  Shown  for  comparison  (as  a  dashed  line) 
is  the  contribution  from  the  first  internal  tide  mode  at  Ml  (Fig.  5-8). 
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Figure  5-12:  RMS  of  the  travel  time  perturbations,  using  the  first  three  IW  modes, 
weighted  by  0.58,  0.26,  and  0.16,  and  from  /  to  N.  The  dashed  line  is  the  rms  of 
the  travel  time  perturbation  for  the  M2  component.  The  dotted  line  is  rms  of  the 
travel  time  perturbation  for  the  across  shelf  case.  For  segment  #2  the  dotted  line  lies 
directly  over  the  dashed  line. 
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slightly  less  than  by  waves  going  along  shelf  (8  =  90),  since  for  the  long  wavelength 
of  the  first  A/2  mode  (260km),  the  sine2  function  has  a  value  of  0.967  as  opposed  to 
1.0  along  shelf.  (Thus  internal  wave  scattering  is  roughly  independent  of  6\)  For  the 
internal  waves,  looking  at  10  cpd  as  a  representative  frequency,  the  across  shelf  sine2 
function  is  0.0001  (compared  to  1.0  for  along  shelf)  due  to  the  higher  frequency  and 
smaller  wavelength.  This  compounds  the  order  of  magnitude  difference  in  spectral 
amplitudes. 

5.3      Movement  of  the  Front 

It  was  pointed  out  earlier  in  Section  3.1  that  the  BSPF  can  move  up  to  10km  over  a 
tidal  cycle.  To  first  order,  this  change  can  be  approximated  by  taking  the  differences 
of  the  average  sound  speed  between  each  segment,  acting  over  one  kilometer.  The 
result  is  a  travel  time  perturbation  of  4.9  msecs/km,  or  up  to  49  msecs  over  a  tidal 
cycle.  As  the  acoustic  array  is  not  expected  to  be  able  to  filter  out  any  modes  higher 
than  mode  #5,  the  modal  travel  time  perturbations  calculated  in  Sections  5.1  and 
5.2  are  not  expected  to  be  a  significant  source  of  noise  for  modal  inverses  for  frontal 
structure.  Internal  wave  noise  could  be  significant,  however,  for  the  ray  (higher  mode) 
arrivals,  a  subject  needing  further  investigation.  Also,  the  phase  relationship  between 
the  perturbation  in  travel  time  caused  by  the  internal  tide  and  internal  wave  fields 
and  the  movement  of  the  front  is  unknown,  because  of  the  complex  relationships 
involved  in  the  bathymetry.  Finally,  due  to  the  fact  that  one  does  have  a  fair  degree 
of  horizontal  resolution  in  the  cross  front  acoustic  maps  made  tomographically,  we 
may  not  be  faced  with  a  significant  ambiguity  between  the  vertical  movement  of 
the  internal  waves  and  the  horizontal  movement  of  the  front.  That  is,  the  internal 
tide/internal  wave  "noise"  may  not  significantly  degrade  the  maps  being  made  of  the 
front. 
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Chapter  6 


Conclusions  and  Future  Work 


As  stated  in  the  beginning,  our  goal  in  this  thesis  was  to  examine  using  acoustic  to- 
mographic techniques  to  study  internal  waves  and  internal  tides.  Towards  that  end, 
we  formulated  an  expression  for  the  travel  time  perturbation  of  an  acoustic  wave  by 
the  internal  tide  and  internal  wave  fields.  The  perturbation  is  maximum  when  the 
acoustic  wave  fields  and  internal  wave  fields  are  propagating  at  right  angles  to  each 
other.  Using  temperature  records  at  three  deployed  moorings,  along  with  the  tem- 
perature profiles  from  CTD  casts,  particle  displacement  records  were  obtained.  We 
then  determined  the  internal  tide  and  internal  wave  spectra  using  Fourier  analysis 
techniques.  The  shape  of  the  particle  displacement  amplitudes  at  the  SW  mooring 
allowed  us  to  estimate  the  amplitude  weights  of  the  first  three  internal  wave  modes. 
Knowing  the  amplitude  spectra  we  then  predicted  the  acoustic  travel  time  pertur- 
bations. We  found  that  for  internal  wave  modes  of  equal  amplitude  that  the  third 
internal  wave  mode  has  the  most  effect  on  the  acoustic  modes.  For  the  Ml  tide,  the 
first  internal  wave  mode  caused  most  of  the  travel  time  perturbation  and  variance. 
When  the  internal  wave  band  was  added,  little  change  from  the  M2  result  was  seen 
for  the  first  five  acoustic  modes.  Since  the  acoustic  VLA  is  not  expected  to  be  able 
to  resolve  modes  higher  than  the  fifth  acoustic  mode,  for  this  application  the  use  of 
the  first  internal  wave  mode  at  the  M2  tidal  frequency  contains  the  bulk  of  the  travel 
time  perturbation  and  variance.  We  next  compared  the  magnitudes  of  the  travel  time 
perturbation  to  the  estimated  change  in  the  front  due  its  movement.  We  concluded 
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that  since  the  front  could  move  upwards  of  10  km  during  a  tidal  cycle,  which  would 
result  in  a  travel  time  perturbation  of  up  to  49  msecs,  that  the  travel  time  perturba- 
tions due  to  the  internal  tide  and  internal  wave  fields  would  be  a  small,  second  order 
effect. 

Whereas  in  this  thesis  we  have  used  Fourier  analysis  techniques  to  extract  the 
internal  tide  and  internal  wave  amplitudes,  there  exist  several  other  methods  which 
could  be  employed  with  varying  degrees  of  success.  One  method  is  the  Lagrangian 
contour  following  method  used  by  Wang  et  al.  [1991]  that  traces  the  minimum  tem- 
perature through  vertically  interpolated  temperature  records.  Another  method  is  to 
determine  the  vertical  velocities  using  the  principle  of  the  conservation  of  temperature 
which  Holloway  [1984]  gives  as 

dT/dt  ,      x 

dT/dz  [      ' 

Both  methods  are  more  suited  to  situations  where  there  are  several  temperature 
records  in  the  vertical  and  where  a  minimal  amount  of  smoothing  needs  to  be  done 
to  remove  the  "spiking"  due  to  eddies  and  other  oceanographic  processes.  There  also 
exist  more  sophisticated  methods  of  tidal  harmonic  analysis  which  can  be  explored 
to  extract  the  amplitudes  such  as  those  used  by  Foreman  [1977]  (cited  in  Holloway, 
1984).  Finally,  a  last  method  to  estimate  the  height  of  the  M2  internal  tide  would 
be  to  calculate  the  movement  of  the  thermocline  using  the  anchor  station  CTD  casts 
at  the  SW  array. 

The  analysis  was  hindered  by  our  ignorance  of  the  direction  of  propagation  of  the 
internal  wave  field.  Levine  and  Richman  [1989]  have  devised  a  method,  whereby  using 
a  single  mooring  of  current  meters  and  thermistors,  the  direction  could  be  resolved. 
Perhaps  such  an  analysis  of  both  the  current  meter  and  the  temperature  information 
may  result  in  resolving  the  propagation  direction.  This  information  would  allow  us 
to  check  our  calculations  versus  the  measured  data  with  much  less  uncertainty. 

Much  more  work  remains  to  be  done  in  fully  understanding  the  nature  of  the 
internal  tide  and  internal  wave  fields.  The  internal  wave  spectrum  could  be  compared 
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to  the  Garrett-Munk  internal  wave  model.  The  travel  time  perturbation  covariance 
was  computed  for  the  0  =  0  and  6  =  90  cases,  but  could  also  be  computed  for  the 
isotropic  case.  Finally,  future  work  in  the  Barents  Sea  could  be  done  to  resolve  the 
phase  relationship  of  the  the  displacement  of  the  BSPF  to  the  barotropic  tide  and 
to  the  internal  tide,  thereby  allowing  a  more  accurate  determination  of  the  acoustic 
signal  travel  time  perturbation. 
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